Medium effects on the quarkonia states above critical temperature 
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We present the quarkonia correlators for charmonium and bottomonium systems in the pseu- 
doscalar, vector and scalar channels. For the description of quark-antiquark interaction we adopt 
the temperature dependant colour screening potential of the power law form. The spectroscopic 
parameters defined from the model are employed in the spectral functions to compute the quarko- 
nia correlators. We find considerable medium modifications to the effective masses of the quarkonia 
as well as in the behaviour of the respective radial wave functions. These modifications are then 
reflected in the computed correlators. The general behaviour of correlators in the vector and scalar 
channel are in accordance with the latest lattice results while their behaviour in the pseudoscalar 
channels are found to be different. 



I. INTRODUCTION 

The behaviour of heavy quarkonium states (charmo- 
nium and bottomonium) in a hot and dense medium have 
attracted considerable experimental and theoretical in- 
terest. Heavy quarkonia play an important role in study- 
ing hot and dense strongly interacting matter. Study of 
the properties of heavy quarkonia above deconfinement 
temperature is of extreme interest for current Relativistic 
heavy ion collider (RHIC) experiments [Ji]. The top- 
ics of interest in these studies include survival probabili- 
ties as bound state at some temperature in quark gluon 
plasma and in medium transport properties of heavy and 
light quarks. Because of larger quark mass the heavy 
quarkonia can survive and remain in bound state above 
the deconfinement temperature T c . The binding energy 
analysis has become important for such studies based on 
either potential models Q or by Lattice methods 
The lattice results of quarkonia correlator studies pre- 
dicts existence of 15 charmonium states upto 1.6T C and 
of IP states (xco and Xci) states upto 1.1T C . The present 
results on bottomonium study predicts existance of IS 
states (rjb and T) upto temperatures 2.3T C and for IP 
states (xb) upto 1.15T C . Which are in contradiction with 
the earlier potential model calculations which predict dis- 
solution of charmonium states below 1.1 T c [7|, |8j. With 
the large number of excited charmonia and bottomo- 
nia states known our earlier understanding based on po- 
tential models have been reviewed recently [3| and sys- 
tematic variations in the confinement strength has been 
introduced to explain the observed quarkonia spectra. 
Such a modification to the string tension of the confin- 
ing part of the potential then be understood in terms of 
the medium effects. Based on such an attempt, here we 
extend our earlier works on the spectroscopy and decay 
properties of quarkonia using coulomb plus power law po- 
tential (CPP V ) model to a systematic study based on the 



quarkonia correlators for charm and bottom quark sys- 
tems. To incorporate the medium effects on the binding 
energy of quarkonia we generalize the usual coulomb plus 
power law (CPP V ) form of the potential with a medium 
dependant exponential screening factor which reduces to 
the coulomb plus power law form in the absence of the 
medium effect. In section |n] we discuss the method to 
compute spectral functions and correlators. The spin av- 
erage mass and the wave function dependant decay con- 
stants are derived in section Hill Finally, we discuss and 
summarize our results for quarkonia correlators against 
different choices of power exponent v in section |V1 



II. EUCLIDEAN CORRELATORS AND 
SPECTRAL FUNCTIONS 

The temperature dependence of the meson correla- 
tors can provide information about the fate of quarko- 
nia states above deconfinement. The imaginary time Eu- 
clidean correlation functions of meson currents G(r, T) 
are reliably calculated on the lattice 



G(t,T) = J du;a{uj,T)K(T,uj,T) 



(1) 



Where, er(w,T) is the zero temperature spectral function 
and K(t,ui,T) is the kernel of integration and can be 
written as, 



K(T,L0,t) 



COsh(uj(T — 



(2) 



In lattice QCD calculation the spectral function <t(uj) can 
be extracted out from the information of correlators us- 
ing the maximum entropy method [9]. While, the poten- 
tial models employ the spectral function to extract out 
the quarkonia correlators. Following ref [l(| the spectral 
function can be written as, 
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The first term arises from the pole contributions from 
the bound states, and the second term is the perturbativc 
continuum above some threshold sq. Here, we consider 
the form of /(w, sq) motivated by leading order pertur- 
bative calculations with massive quarks as Q, 



/(w,So) 



a H 



(4) 



The calculated coefficients of an and bn at leading or- 
der by [ll| are shown in Table [fl The value of threshold 
energy so is somewhat arbitrary for single flavour heavy 
quark only. Here, Sq is chosen such that no resonance 
above it is possible. The remaining parameters Mi and 
wave function dependant decay constants Fi are deduced 
from the potential model adopted for the present study. 
To see the temperature effect on the spectral function and 
to compare with the lattice QCD results one usually com- 
putes the ratio of this correlators to the reconstructred 
one G(r, T)/G 

r i CO!/ 

(r, T), where G r (r, T) is given by, 



>- 1.60 




1.0 1.2 



FIG. 1. V(r — oo,T) for different temperature at different 
potential exponent 



G™(r,T)= / d W( r( W ,r = 0)^(T, W ,r) (5) 
Jo 

Here, G recon (T, T) corresponds to the spectral function 
at the zero temperature. The temperature dependance 
in the G recon (t, T) comes only from the kernel of the 
integration. 

TABLE I. The coefficients (an, bn) in different mesonic chan- 
nel m 



System 


a H 


b H 


scalar 


-1 


1 


pseudoscalar 


1 





vector 


2 


1 


axial- vector 


-2 


3 



III. EXTRACTION OF SPIN AVERAGE MASS 
(Mi) AND DECAY CONSTANTS (Fi) 

The spin average mass and decay constants appeared 
in the expression of the spectral function (Eqn. [3|) are 
deduced using an appropriate model description of the 
hadronic state. Here, for the description of the quarko- 
nia states we consider temperature dependant screened 
potential of the form, 

V(r,T) = + _£_(! - e-AfW) ( 6 ) 

Here, a and a are the coupling constant and string 
tension respectively. The exponent, v is the variation 
that we have introduced here to account for the non- 
linear dependance on the interquark distance of the po- 
tential. Different choices of the exponent, v describes 



different strength of the interquark interaction. The 
screening mass parameter w(T) is taken as, fi(T) — 
0.24 + 0.31 (T/T c - 1) GeV Q with critical temperature 
T c = 0.270 GeV. Here, we have chosen A = 1 GeV"- 1 . 
Here, the screening potential attains finite value at in- 
finite separation i.e. V(r — oo,T). Hence, each quark 
has an additional thermal energy of Vqq/2 in the bound 
state. The minimum energy above which the quark- 
antiquark pair can freely propagate is 2m + Voq. Hence, 
the minimum energy threshold can be defined in this 
case as sq — 2m + V^. The potential Voo for differ- 
ent temperature is shown in Fig. [T]for different potential 
exponent. Similarly, the pole mass Wj corresponds to 
mi = m + Voo/2 against T/T c for different potential ex- 
ponent v are shown in Fig. [3 for bottomonium and char- 
monium states. Voo found to be decreasing with increase 
in temperature above T c . Similar observation is observed 
on lattice calculations while calculating quark and gluon 
propagators in coulomb gauge [13 • In the absence of 
the medium effect (at zero temperature) the potential in 
Eqn. [6] reduces to the form 

V(r) = -- + <rr u (7) 
r 

The parameters a, a and the quark masses are fixed 
by fit ting the zero temperature quarkonium spectrum 
as in jf 3l . [lij]. The parameters employed are given by 
a = 0.471, a = 0.192 GeV v+1 , m c = 1.32 GeV and 
m& = 4.746 GeV by fitting the zero temperature quarko- 
nia spectrum. The spin average masses and the wave 
functions are obtained by solving the schroedinger equa- 
tion numerically using the screened potential of Eqn. [6] 
as quark- antiquark interaction potential at a given tem- 
perature. 

In the relativistic quark model, the decay constant can 
be expressed through the meson wave function <& P / V {p) 
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FIG. 2. The pole mass for bottom quark for screened coulomb potential at different potential exponent (a) for bb system (b) 
for cc system 
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FIG. 3. rms deviation in zero temperature mass for quarkonia states taken from [13] (a) for charmonium states (b) for 
bottomonium states 



in the momentum space [If 
d 3 p 



fp/ 




E Q (p) + m Q \ ( E Q (p)+m 



(2tt) 3 \/ V 2E Q (p) 
[E Q (p) + m Q ][E Q (p) + m Q ] 



2E Q (p) 
®p/v(p) 



with Ap = —1 and Ay- = —1/3. In the nonrelativis- 

2 

tic limit ^2 << 1.0, this expression reduces to the well 
known relation between fp/v and the ground state wave 
function at the origin R P / V (0) through the Van-Royen- 
Weisskopf formula 16]. Though most of the models 
predict the meson mass spectrum successfully, there ex- 
ist wide range of predictions of their decay constants. 
For example, the ratio was predicted to be > 1 



in most of the nonrelativistic cases, as mp < my and 
_thcir wave function at the origin assumed to be the same 
i V?p(0) ~ Ry (0) [I3]- The ratio computed in the relativis- 
ytic 



tic models [18| have predicted 4 s - < 1, particularly in the 



fv 

^ sector, but 4 s - > 1 in the heavy-light flavour sector. 
(&£ ■ ■ * v 

v This disparity on the predictions of the decay constants 
play decisive role in the decay properties of these mesons. 
The value of the radial wave function (Rp) for and 
(Rv) for 1 states would be different due to their spin 
dependent hyperfine interaction. The spin hyperfine in- 
teraction of the heavy flavour mesons are small and this 
can cause a small shift in the value of the wave function 
at the origin. Though, many models neglect this differ- 
ence between (Rp) and (Rv), we consider this correction 
by making an ansatz that the i?p/y(0) are related to the 
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FIG. 4. Ratio of the mass at given temperature to the zero temperature mass for different quarkonium states at different 
potential exponent v (a) for charmonium (b) for bottomonium. 




value of the radial wave function at the origin, i? n (0) ac- 
cording to the same way their masses are related. Thus, 
by considering 



M nP/v = M„ 



n,CW 



1 + (SF) Pn 



(Vss)r 



M, 



n,CW 



(9) 



and following the fact that any c-number, a, commutes 
with the Hamiltonian, i.e. aHty = H(a^), we express 

EEi, 



RnP/v(0) = Rn{0) 



1 + (5F) P/1 



(M nV - M nI 



cw 



(10) 



Here (SF)p = — | and (SF)y — j are the spin factor 
corresponding to the pseudoscalar (J = 0) spin coupling 



and vector (J = 1) spin coupling respectively. M nt cw 
and i? n (0) are spin average mass and the normalized spin 
independent wave function at the origin of the meson 
state respectively. It can easily be seen that this expres- 
sion given by Eqn [10] is consistent with the relation 



R(0) = 



3R V (0) + Rp(0) 



(11) 



given by [19|,l2Q| for nS states. By incorporating first or- 
der QCD correction to the Van Royen-Weiskopff formula, 
the decay constant is computed as [2l|, [22j , 
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p/v 



-C\a s ) 



(12) 
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FIG. 6. G/Grecon for pseudoscalar charmonium (j] c ) state (a) for v — 1.0 (b) for v — 1.1 
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FIG. 7. G/Grecon for vector charmonium (J/tp) state (a) for ^ = 1.0 (b) for f = 1.1 



where, the first order QCD correction, C(a s ) is expressed 
for the QQ system as 



C(a s ) = 1-^-8 

IT 



(13) 



Here S v = f |2l|, [22] and 5 P = 2 [aM 



IV. TEMPERATURE AND THE EXPONENT (y) 
DEPENDANCE ON THE CORRELATORS 



The temperature dependant energy eigen values 
(i?i)and the wave functions are obtained by solving the 
schrodinger equation numerically [24| . The spin average 
mass (M sa ) is given by M sa = 2irii + Ei. Once M sa is 
obtained, then the mass of the quarkonia state is com- 
puted using Eqn. HI For computing the mass difference 



between different spin degenerate mesonic states, we con- 
sider the spin dependent part of the usual one gluon ex- 
change potential (OGEP) given by [2514291 ]. Accordingly, 
the spin-dependent part, Vsd(t) contains three types of 
interaction terms, such as the spin-spin, the spin-orbit 
and the tensor part as 



Vsd(t) = V S s(r) 
V T (r) 



5(5- 
5(5 4 



1 



1) - 



V LS {r) (L-S) 



3 
" 2 

3(5-f)(5-f) 



(14) 



The spin-orbit term containing Vi,s(r) and the tensor 
term containing Vr{r) describe the fine structure of the 
meson states, while the spin-spin term containing Vss(r) 
proportional to 2(s^ • s^) = 5(5 + 1) — | gives the spin 
singlet-triplet hyperfine splitting. The coefficient of these 
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FIG. 8. G/Grecon for scalar charmonium state (Xco) (a) for v — 1.0 (b) for v = 1.1 
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FIG. 9. G/Grecon for pseudoscalar bottonium state (??{,) (a) for v — 0.7 (b) for v — 1.0 



TABLE II. Spin averagem mass (M sa ) in GeV for charmonium and bottomonium states for different choices of potential 
exponent v 

Charmonium state Bottomonium state 



IS IP IS IP 



T/T c 


v =0.5 


1.0 


1.1 


1.5 


0.5 


1.0 


1.1 


1.5 


0.5 


0.7 


1.0 


1.5 


0.5 


0.7 1.0 


1.5 


1.1 


2.930 


3.040 


3.060 


3.129 


3.063 


3.274 


3.306 


3.370 


9.515 


9.519 


9.525 


9.535 


9.787 9.831 9.891 


9.976 


1.2 


2.929 


3.034 


3.053 


3.115 


3.053 


3.242 


3.268 


3.305 


9.523 


9.527 9.533 


9.542 


9.787 


9.829 9.887 9.967 


1.3 


2.928 


3.028 


3.045 


3.100 


3.042 


3.211 


3.229 


3.254 


9.530 


9.534 


9.540 


9.549 


9.786 


9.827 9.883 


9.958 


1.4 


2.926 


3.021 


3.036 


3.085 


3.033 


3.181 


3.195 


3.212 


9.537 


9.541 


9.547 9.555 


9.784 9.824 9.878 


9.947 


1.5 


2.923 


3.013 


3.027 3.070 


3.024 


3.154 


3.165 


3.179 


9.544 


9.548 


9.554 


9.562 


9.782 


9.821 9.872 


9.935 


1.6 


2.921 


3.005 


3.018 


3.054 


3.018 


3.132 


3.141 


3.153 


9.550 


9.554 


9.560 


9.568 


9.780 


9.817 9.865 


9.922 


1.7 


2.918 


2.997 


3.008 


3.039 


3.013 


3.113 


3.121 


3.131 


9.556 


9.560 


9.566 


9.573 


9.776 


9.813 9.858 


9.907 


1.8 


2.915 


2.989 


2.999 


3.025 


3.009 


3.099 


3.105 


3.114 


9.562 


9.566 


9.572 


9.579 


9.773 


9.808 9.850 


9.892 


1.9 


2.912 


2.980 


2.989 


3.012 


3.007 


3.087 3.093 


3.101 


9.568 


9.572 


9.577 


9.584 


9.769 


9.803 9.842 


9.876 


2.0 


2.910 


2.973 


2.981 


3.000 


3.007 


3.079 


3.084 


3.091 


9.573 


9.577 


9.582 


9.589 


9.766 


9.797 9.834 


9.862 
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FIG. 10. G/Grecon for vector bottonium state (T) v — 0.7 (a) for v = 0.7 (b) for v = 1.0 
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FIG. 11. G/Grecon for scalar bottonium state (xf>o) (a) for v = 0.7 (b) for v = 1.0 



TABLE III. Masses of charmonium states in GeV for different choices of potential exponent v 



T/T c v=0.5 1.0 1.1 1.5 0.5 1.0 1.1 1.5 0.5 1.0 1.1 



1.5 



1.1 2.886 2.960 2.973 3.020 2.945 3.067 3.089 3.166 

1.2 2.888 2.960 2.973 3.017 2.943 3.059 3.079 3.148 

1.3 2.889 2.959 2.972 3.014 2.940 3.051 3.069 3.129 

1.4 2.890 2.958 2.970 3.010 2.938 3.041 3.058 3.110 

1.5 2.890 2.957 2.968 3.006 2.934 3.032 3.047 3.091 

1.6 2.890 2.955 2.965 3.000 2.931 3.022 3.035 3.072 

1.7 2.890 2.952 2.962 2.993 2.927 3.012 3.024 3.054 

1.8 2.889 2.948 2.957 2.985 2.924 3.002 3.012 3.038 

1.9 2.887 2.944 2.952 2.977 2.921 2.993 3.002 3.023 
2.0 2.886 2.940 2.947 2.969 2.918 2.984 2.992 3.010 



4.312 4.412 4.469 5.482 
4.323 4.464 4.549 5.930 
4.332 4.523 4.639 6.197 
4.337 4.578 4.716 6.322 
4.337 4.619 4.770 6.350 
4.334 4.645 4.801 6.318 
4.327 4.656 4.814 6.249 
4.319 4.657 4.812 6.160 
4.309 4.650 4.800 6.059 
4.297 4.637 4.781 5.955 
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TABLE IV. Masses of bottomonium states in GeV for different choices of potential exponent v 

Vb T XM 

T/T c v =0.5 0.7 1.0 1.5 0.5 0.7 1.0 1.5 0.5 0.7 1.0 1.5 



1. 


.1 


9 


.464 


9 


.464 


9 


.464 


9 


.465 


9 


.532 


9 


.538 


9 


.546 


9 


.558 


II 


.359 


11 


.359 


11 


.360 


11 


.361 


I 


.2 


9 


.472 


9 


.472 


9 


.472 


9 


.474 


9 


.540 


9 


.545 


9 


.553 


9 


.565 


11 


.359 


11 


.359 


11 


.360 


11 


.361 


I 


.3 


9 


.481 


9 


.480 


9 


.481 


9 


.482 


9 


.547 


9 


.552 


9 


.560 


9 


.571 


11 


.359 


11 


.359 


11 


.360 


11 


.361 


I 


.4 


9 


.489 


9 


.488 


9 


.489 


9 


.490 


9 


.554 


9 


.559 


9 


.567 


9 


.577 


11 


.359 


11 


.360 


11 


.360 


11 


.362 


I 


.5 


9 


.496 


9 


.496 


9 


.497 


9 


.498 


9 


.560 


9 


.565 


9 


.573 


9 


.583 


11 


.360 


11 


.360 


11 


.360 


11 


.363 


I 


.6 


9 


.504 


9 


.504 


9 


.504 


9 


.505 


9 


.566 


9 


.571 


9 


.579 


9 


.589 


11 


.360 


11 


.360 


11 


.361 


11 


.363 


I 


.7 


9 


.511 


9 


.511 


9 


.511 


9 


.512 


9 


.572 


9 


.577 


9 


.584 


9 


.594 


11 


.360 


11 


.360 


11 


.361 


11 


.365 


I 


.8 


9 


.517 


9 


.518 


9 


.518 


9 


.519 


9 


.577 


9 


.582 


9 


.589 


9 


.599 


11 


.360 


11 


.360 


11 


.361 


11 


.366 


I 


.9 


9 


.524 


9 


.524 


9 


.525 


9 


.526 


9 


.582 


9 


.587 


9 


.594 


9 


.603 


11 


.360 


11 


.361 


11 


.362 


11 


.367 


2. 


.0 


9 


.530 


9 


.531 


9 


.531 


9 


.532 


9 


.587 


9 


.592 


9 


.599 


9 


.608 


11 


.360 


11 


.361 


11 


.362 


11 


.369 



TABLE V. I th derivative of radial wave function at zero separation for charmonium states in GeV l+3 ^ 2 for different choices of 
potential exponent v 

Vc J_hP_ xco 

r/r e i/=o.5 l.o l.i 1.5 0.5 l.o 1.1 1.5 0.5 1.0 1.1 1.5 



1 


.1 


0. 


239 





.423 


0. 


456 





.565 


0. 


249 





.454 





.492 





.620 





.051 





.102 





.106 





.047 


1 


.2 


0. 


224 





.392 


0. 


421 





.509 


0. 


232 





.419 





.452 





.554 





.046 


0. 


.083 





.080 





.033 


1 


.3 


0. 


208 





.362 


0. 


.387 





.453 


0. 


216 





.385 





.412 





.488 





.042 





.066 





.060 





.029 


1 


.4 


0. 


193 





.331 


0. 


352 





396 


0. 


200 





.350 





.373 





.423 





.039 





.054 





.048 





.028 


1 


.5 


0. 


179 





.301 


0. 


.317 


0. 


.341 


0. 


184 


0. 


.317 


0. 


.334 


0. 


.361 


0. 


.037 


0. 


.047 





.042 





.028 


1 


.6 


0. 


165 





.271 


0. 


.283 





.290 


0. 


170 





.284 





.296 





.304 





.037 





.043 





.040 





.029 


1 


.7 


0. 


154 





.243 


0. 


.251 





.247 


0. 


158 





.253 





.262 





.257 





.037 





.042 





.039 





.030 


1 


.8 


0. 


144 





.218 


0. 


223 


0. 


.212 


0. 


147 


0. 


.226 


0. 


.231 


0. 


.220 


0. 


.038 


0. 


.041 





.039 





.032 


1 


.9 


0. 


136 





.197 


0. 


200 





.187 


0. 


139 





.203 





.206 





.193 





.040 





.042 





.041 





.035 


2 


.0 


0. 


130 





.180 


0. 


181 





.169 


0. 


133 





.185 





.187 





.174 





.042 





.044 





.042 





.038 



spin-dependent terms of Eqn|14] can be written in terms 
of the vector (W) or coulomb part and scalar (Vs) or 
confinement par t of the static potential, V(r) described 
in Eqn. [Bias [27| 

„„ M __L_ (15) 

2 m\Tn,2 r \ ar ar J 

vm 1 fo<PVv ldv v\ n ,, 

V T (r) = 3— — (16) 

m\rri2 \ dr* r ar J 

V SS (r) = ^^S^(f) (17) 

The spin average masses (M sa ) for both charmonium 
and bottomonium states are tabulated in Table 111 The 
quarkonia masses for both charmonium and bottomo- 
nium states are shown in Table IIIII and IIVI respectively. 
Similarly the £ th derivative of the radial wave function at 
zero separation for charmonium and bottomonium states 
are shown in Table IVl and I VII respectively. The choices of 
the potential exponent v = 1.1 for charmonium state and 



v = 0.7 for bottomonium state correspond to the mini- 
mum of the standard deviation in the mass (SD m in) hi 
zero temperature potential as observed in our previous 
study |13| as shown in Fig. [3j The temperature de- 
pendant masses of the quarkonia states in terms of their 
zero temperature masses (Mj(T)/Mj(0)) are computed 
for T > T c and for different choices of the potential ex- 
ponent, v are shown in Fig. |4j similarly the ratio for 
the wavefunction is shown in Fig. [5] The ratio in mass 
for the P-state is found to be larger compare to the S- 
state quarkonia states. The ratio in the wave function 
decreases with increasing tempearature for all quarko- 
nium states. Finally, the decay constants are computed 
with help of Eqn. Ql 

Our computed Spin average masses (M sa ) and the decay 
constants (f p / v ) corresponding to the choises of v = 1.1 
in the charmonia and v = 0.7 in the case of bottomonia 
are then employed to predict the spectral functions given 
in Eqn. [3l We have also computed spectral function for 
the choices v = 1.0 in both the cases for the purpose 
of comparison with cornel like interactions considered by 
others [3j, |3fJ . Further the euclidean correlators are com- 
puted using Eqn. [1] The present ratio G/G recon are 
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TABLE VI. I th derivative of radial wave function at zero separation for bottomonium states in GeV e+3 ^ 2 for different choices 
of potential exponent v 

Vb T XbO 



T/T c v=0.5 0.7 1.0 1.5 0.5 0.7 1.0 1.5 0.5 0.7 1.0 1.5 

1.1 6.810 7.385 8.147 9.182 6.909 7.502 8.289 9.363 0.527 0.671 0.874 1.168 

1.2 6.695 7.258 8.001 9.005 6.791 7.370 8.138 9.179 0.498 0.635 0.826 1.089 

1.3 6.575 7.125 7.852 8.826 6.668 7.234 7.983 8.993 0.468 0.598 0.776 1.006 

1.4 6.451 6.990 7.698 8.642 6.540 7.095 7.825 8.802 0.437 0.561 0.725 0.917 

1.5 6.322 6.851 7.542 8.455 6.407 6.951 7.663 8.608 0.407 0.524 0.672 0.822 

1.6 6.189 6.707 7.381 8.266 6.271 6.803 7.497 8.412 0.376 0.485 0.618 0.715 

1.7 6.054 6.560 7.217 8.072 6.132 6.652 7.328 8.211 0.347 0.447 0.563 0.597 

1.8 5.915 6.409 7.051 7.876 5.989 6.497 7.157 8.008 0.319 0.410 0.506 0.483 

1.9 5.771 6.257 6.881 7.676 5.841 6.340 6.981 7.802 0.294 0.374 0.450 0.392 
2.0 5.626 6.100 6.707 7.475 5.693 6.179 6.803 7.594 0.273 0.342 0.398 0.331 



0.3 



0.1 



■ 








■ Other than 1S bb 
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• 1Sbb 
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■ 
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■ 
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1.0 1.2 1.4 1.6 1.8 2.0 



T/T c 

FIG. 12. T max vs. T/T c for all quarkonia states 



drawn in Figs. EJ H El [TO] and HU for the euclidean 
time t in between r = 0tor = l/T. The plots of these 
figures correspond to different choice of T/T c < 2.0. 



V. RESULTS AND DISCUSSION 

As from Fig. H Mi(T)/Mi(0) for J/ip and rj c case 
found to decrease continuously with T/T c for all the 
range of T/T c . The behaviour is same for both the 
choices of v = 1.0 and v = 1.1 except the magni- 
tude of Mi(T)/Mi(0) increases slightly with v. While 
Mi(T)/Mi(0) for Xco found to increase uptoT/T c = 1.7 
then, start slightly decreasing. However, in the case of bb 
system Mi{T)/Mi(Q) for all T and Xbo states found to 
increases with T/T c for all T in the range of T/T c < 2.0. 
So, the medium effect in the r\ c and J/ip are opposite to 
that of ijb and T states. 

From, Fig. MRi(T)/Ri(0) for rj c and J/ip states found 



to decrease monotonically with increase in T/T c , while 
Rf\T)/Rf\Q) for Xco state found to decrease upto 
T/T c = 1.7 then it is increasing slightly. R t {T)/Ri{0) 
for r/b and T states found to decrease monotonically for 
all the range of T/T c . Rf ] {T)/Rf } (0) for xw state also 
found to decrease much faster than rjb and T states. 
These observations are same for both v = 1.0 and v = 0.7 
except their magnitudes are higher. 
The correlators at the different T/T c obtained are found 
to have a similar behaviour with increase in r. In the case 
of charmonia state, correlators attain maxima at differ- 
ent values of t and it shilfts towards lower r as T/T c 
increases. While in the case of r\b and T bottomonia the 
maxima of the correlators are found to shift slightly to- 
wards higher r. Whereas, the behaviour in the case of 
Xbo belongs to that of charmonia states. To make this 
boservation more clear we plot the values of r at the 
maximum correlator (tau ma x) against T/T c in Fig. [12] 
for both the case of bb and cc. It is very striking to see 
that values of r at which G/G recon is maximum (T max ) 
is same for rj Cl J ftp, Xco an d Xbo states with a chosen 
T/T c , and T max decreases with T/T c . This behaviour is 
shown graphically in Fig. [T3]Here we also shown the be- 
haviour of T max vs T/T c in the case of 775 and T system. 
It seems that both the curves leading towards a common 
saturated value beyond T/T c > 2.0. Our results do not 
show any major deviation with the potential exponent v. 
However, from Fig. [g](b), we find G/G recon &t T/T c = 1.1 
of x c o decreases from 1.0 as r increases while similar case 
in Fig. [8] (a) for v = 1.0 seems to increase from 1.0 as r 
increases. From Fig. [TU] and Fig. |H] we also notice that 
G/G recon in rjb and T decreases below 1.0 relatively at 
lower r values compared to all other states studied here. 
The r dependance of the correlators are found to be sen- 
sitive to the choices of the continuum threshold so- 
As seen from Fig. [HI our obtained G/G recon shows agree- 
ment with the early potential model Q but differs from 
the lattice results Q where G/G recon remains to unity 
upto t w 0.05 fm then it starts gradually falling. Our 
results shows same trend to QCD sum rule results [3(| at 
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T = 1.1T C . 

Fig. [7] represents G/G recon for vector charmonium states 
(J/ip). For smaller values of r our results agrees with the 
lattice @ and also with the QCD sum rules [3(| ■ Though 
the present G/G recon in the case of rj c increases from 1.0, 
the results based on fine isotropic lattice study shows de- 
creases from 1.0 for all temepratures. However the gen- 
eral behaviour of all other states are in accordance with 
the lattice results As, the r] c and J/ip corresponds 

to same IS state, one do not expect their correlators to 
behave differently. 

Fig. [5]represents G/G recon for scalar channel of IP char- 
monium state (xco)- Here, also our results agrees with 
lattice @ and QCD sum rules (30j . 

Similarly, Figs. [9| ITOl and ITT1 represent our results for IS 
pseudoscalar (rjb), IS vector (T) and IP scalar (xbo) bot- 
tomonium states respectively. Our results agrees with 
earlier potential model Q but contradicts with lattice 
results 5]. The lattice results for rjb case shows unity 
behaviour for temperature upto 2.30T C . The weaker de- 
pendance of the correlators on temperature was assumed 
as the melting of the states, but melting of state can- 
not be exctracted from the correlators [32|]. Similar to 
the charmonium case, our results for rjb and T shows 
the same behaviour as expected (because they belong to 
same IS state of bottomonium system). 
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